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Abstract 

Combining a lattice path integral formulation for thermodynamics with the solution 
of the quantum inverse scattering problem for local spin operators, we derive a multiple 
integral representation for the time-dependent longitudinal correlation function of the 
spin-1/2 Heisenberg XXZ chain at finite temperature and in an external magnetic field. 
Our formula reproduces the previous results in the following three limits: the static, the 
zero-temperature and the XY limits. 

PACS numbers: 05.30.-d, 75.10.Pq, 02.30.1k 



1 Introduction 

One of the most challenging problems in quantum many-body systems, especially in low- 
dimensions, is related to the evaluation of the correlation functions. Though many field 
theoretical schemes or numerical techniques have been developed and achieved remarkable 
success, the exact evaluation is in general still very difficult even in models which are exactly 
solvable by the Bethe ansatz methods. 

The spin-1/2 Heisenberg XXZ chain is one of the simplest but non-trivial solvable mod- 
els, and has served as a testing ground for many theoretical approaches. Concerning the 
correlation function of this model, Jimbo et al. [1,2] derived a multiple integral representa- 
tion of arbitrary correlators in the off-critical XXZ antiferromagnet at zero temperature and 
zero magnetic field. Their method utilizing the representation theory of the quantum affine 
algebra U q {s\2) has been extended to the XXX [3,4], the massless XXZ [5] antiferromagnets. 

Alternative approaches combining the algebraic Bethe ansatz and the solution of the 
quantum inverse scattering problem for local spin operators were proposed by Kitanine et al. 
in [6] , and were applied to derive a multiple integral representation for both the critical and 
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off-critical regimes of the XXZ model in finite magnetic field. Moreover they have obtained a 
new multiple integral representation which is more appropriate for the two-point correlators 
[7, 8] , and have generalized further to the time-dependent longitudinal correlation function 
[9] (see [10] for a recent review). On the other hand, a finite temperature generalization 
was achieved more recently by Gohmann et al. by utilizing the quantum transfer matrix 
formalism for thermodynamics [11] (see also [12-14] for recent progress in this direction). 

Motivated by these seminal works, in this paper, we generalize the result [9, 11] to the time- 
dependent longitudinal correlation function at finite temperature. By combining a lattice 
path integral formulation with the solution of the quantum inverse scattering problem, we 
derive a multiple integral representing a generating function of the time-dependent correlation 
function for the ^-components of the spins. In the zero-temperature, the static and the XY 
limits, our formula reproduces the results in [9, 11, 15]. 

The layout of this paper is as follows. In the subsequent section we introduce the spin- 
1/2 Heisenberg XXZ chain and its classical counterpart. Utilizing a lattice path integral 
formulation together with the solution of the quantum inverse scattering problem for local 
spins, in Section 3, we formulate the generating function characterizing the time-dependent 
longitudinal correlation function at finite temperature. In Section 4, we derive a multiple 
integral representing the generating function by combining the method developed in [7, 9] 
and [11]. The following three limits, namely the static, the zero-temperature and the XY 
limits are considered in Section 5. Section 6 is devoted to a brief conclusion. 



2 Spin-1/2 Heisenberg XXZ chain 



The Hamiltonian of the spin-1/2 Heisenberg XXZ chain defined on a periodic lattice of length 
L is given by 



where 



j'=i 



(2.1) 



(2.2) 



Here aj ,y ' z are the Pauli matrices acting on the two-dimensional quantum space Hj at site 
j. The real parameter A>o is the anisotropy parameter, and J € R fixes the energy scale of 
the model. In this paper, we consider the model at finite temperature T > and in a finite 
magnetic filed h € R. 

It is well-known that a d-dimensional quantum system can be mapped onto a (d + 1)- 
dimensional classical system. The classical counterpart of the XXZ model is so-called the 
six-vertex model whose Boltzmann weights can be described as the elements of the ^-matrix 
&{\) G End(C 2 ®C 2 ), 
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Figure 1: A graphical representation for the ^-matrix (2.3) and the monodromy matrix 

55(A) (2.4). 

Identifying one of the two vector spaces of the above ^-matrix with the quantum space T~Cj , 
we define the monodromy matrix ^ R (A) as 



5f(A) 



HL 



(A) 



(A)^i(A), 



(2.4) 



where ^ 7 (A) acts in the space 7Yj <8> W,- (see also figure 1 for a graphical representation of 



and ^ R (A)). Since the ^-matrix satisfies the Yang-Baxter equation 
^12 (A - m)#13(A - v)&2a(n -") = ^23 (/u - (A - i/)^i 2 (A - /x), 



(2.5) 



the transfer matrix defined by Tr(A) = Tr^ ^ R (A) commutes for different spectral parameters: 
[T R (A), T R (/i)] = 0. Using the following relation connecting the six-vertex model with the 
XXZ model (2.1): 

J^o = 2Jsh(r ? )T- 1 (0)T^(0), A = ch <q, (2.6) 



one can expand the transfer matrix with respect to the spectral parameter A 



Tr(A) = T R (0) 1 + 



A 



-J^o + 0(X 2 



(2.7) 



2Jshr/ 

For later convenience, let us introduce another type of transfer matrix defined by 

T R (A) = Tr; fe(-A) • • • ^(-A)^(-A)] . (2.8) 
Note that Tr(0) and T R (0) are respectively the right- and left-shift operators, and hence 

T R (0) = T- 1 (0). (2.9) 
Using this together with the expansion (2.7), one arrives at a crucial formula 



lim 



T R (A)T R A + 



N 



exp 



A=o 



4J shr] 



-JtU 



(2.10) 



where N € 2N is assumed. 



1 We call the remaining space the auxiliary space and write Ttj 
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3 Time-dependent generating function at finite temperature 
3.1 Quantum transfer matrix 

Our main aim in this paper is to derive a multiple integral representation of the time- 
dependent longitudinal correlation function (erf (0)<r^ +1 (t)) for the XXZ model (2.1) at finite 
temperature. To this end, we firstly consider how the time-dependent correlator is described 
by the transfer matrix formalism. 

Since [S z ,aj] = 0, the time-dependent local spin operator cr z {t) can be written as 

o-|(t) =e im a z j e- i ^ ot . (3.1) 
Using this and noticing that [J^o, S z ] = 0, one immediately sees 

K(0)< +1 (*)> = ^ ^ e _J/ T (3-2) 

Here we have set the Boltzmann constant to unity. Applying the formula (2.10), we insert 
the relations 



e-W= lim e^{T R (A)T R (A- £o )}~ 

No — »oo 



A=0 



e ±im = lim {T R (A)T R (A± £l )} : 

Ni— >oo 



A=0 



so = —, 00 = —^, 

ei = TT, 01 = 4itJshr/, (3.3) 
iVi 



into (3.2). The result can be graphically represented as in figure 2. As proved in [16-18], the 
local spin operator cr| is expressed in terms of the transfer matrix: 

a] = Tl~\0) Trj\a? ^ R (0)]T^(0). (3.4) 

Thus, substituting 

a\ = hm Tr^af ^ R (A - e)]T R (A) e G C (3.5) 

£ *U A — U 

which is obtained by (3.4) and (2.9), one finds (af (0)cr^ +1 (i)) is given by a product of the 
transfer matrix and two local operators located on the boundaries (see figure 3). 

To consider the system at finite temperature, let us introduce the quantum transfer matrix 
T(A) acting in the space ®£+ 2Nl+2 Hr. 

T(A) = Trj ^j(X). (3.6) 

Here ^j(X) is defined by 

x %]vT+2(- A )%vr+T,( A - £l ) • • -^Th+ii-^W+Ej^ ~ £ i) 
x^7vr+2,( A - £ )^]vr+T(- A ) 

x ^-(-A)^i.(A + £!)■•■ %(-A)^ Tj (A + (3.7) 
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Figure 2: A graphical representation for the longitudinal dynamical correlation function 
(<rf : (0)(T^ t+1 (t)). Here = /3q/Nq and e\ = (3\/N\. We assume the lattice is on a torus. The 
correlation function (multiplied by the partition function Tr e~~^l T ) is reproduced by setting 
A = and taking the limit A^o,i — > oo. The partition function Tre~'^/ T is obtained by just 
replacing o\ and cr^ , 1 with the identity operator. 



In figure 3, we also schematically depict the quantum transfer matrix. Thanks to the Yang- 
Baxter equation (2.5), we see the quantum transfer matrix commutes for different spectral 
parameters: [T(A), T(//)] = 0. Thus the dynamical correlation function (3.2) is expressed in 
terms of T(A): 



Tn 



(af(0)a z m+1 (t)) = lim lim- 

7Vo,-/Vi^oo e— >0 



-iL—m—l 



(e)(A-D)(e)T m (e)a. 



m+2 



Tr 



1,...,L 



(3.8) 



where A(X) and D(X) are elements of the monodromy matrix ^(X) represented as a 2 x 2 
matrix 

A(X) B(X) 
C(X) D(X) 



3T{X) 



(3.9) 



in the quantum space. In this definition, T(A) is given by T(A) = A(X) + D(X). 
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Figure 3: A graphical representation of the quantum transfer matrix T(A) = Tr, ^j(X) 
(surrounded by broken lines). In fact, a small parameter e is introduced to avoid that the 
quantum transfer matrix becomes a singular matrix. The partition function Tre~"^/ T is 
obtained by replacing cr^ +1 and c z Ni+2 with the identity operator. 

Let us consider the thermodynamic limit L — > oo. Since the limits L — > oo and Nq,N\ — > 
oo are interchangeable as proved in [19, 20], one can take the limit L — > oo first. In addition, 
we find the leading eigenvalue of the quantum transfer matrix T(0) (written as Ao(0)) is non- 
degenerate and separated from the next-leading eigenvalues by a finite gap even in the Trotter 
limit Nq,Ni — > oo. In the thermodynamic limit L — > oo, therefore, (3.8) is characterized by 
Ao(0) 2 and the corresponding (normalized) eigenstate |^): 

K(O)CiW) = w lim n A o m " 'OOO^P - J D)(e)T m (e)^|v&). (3.10) 

Inserting the relation 

°k+* = [T(si)T-\0)] ^ T-\0)(A - D)(e)T~\e)T(0) [T^T'^-e,)] ^ , (3.11) 
2 Note that A (A) is defined by A (A) = (*|T(A)|*>. 
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which is a 'quantum transfer matrix analogue' of (3.4), we finally obtain 
K(0)< + iW) = M Hm lim(*|(^ - D){e)T m ~\e)Q 

N ,Ni— >oo£— >0 

x (A - D)(e)Q" 1 T" m " 1 (e)|^), (3.12) 

N i 

where Q = T(e)T- 1 (0) [T- 1 (0)T(-e 1 )] 2 . 

3.2 Generating function for dynamical correlation function 

It is convenient to introduce the following operator as in [9]: 

£i+i,m = T l {e)T™-\e)Q R T- m (e)Q-\ (3.13) 
Here T K (X) and Q K are respectively defined by 

T K (X) = A(X) + kD(X), Q k = T K (e)T-\0) [T- 1 (0)T K (-e 1 )] ^ . (3.14) 



Due to the Yang-Baxter equation (2.5), the twisted quantum transfer matrix T K (X) is com- 
mutative as long as the twist angle k is taken the same: [T K (X),T K (fi)] = 0. It follows that 

= , m i im( ^|^ (£)jD(£)r -^i (£)Q 

2 1 / Nq,Ni—>oo e^O 

x £>(e)Q- 1 T- m - 1 (e)|*) 

= N J^J™ \dlm^l + l,m + l ~ £l + 2,m + l - £l + l,m + ^f+2,J*> \ K=V (3-15) 

Because of the translational invariance for the correlation functions, one can set I = without 
loss of generality. Introducing the generating function 

Q K (m,t) = hm lWtfl^Jtt), (3.16) 

iVo,iVl— >oo £^0 ' 

one obtains the longitudinal time-dependent correlation function (3.2): 

W(0Km+iW) = 2^^Q K (m,t)| K=1 +2(^) - 1, (3.17) 

where (cr 2 ) is the magnetization (multiplied by a factor 2) and D m denotes the lattice deriva- 
tive defined by 

D m g{m) = g{m + 1) - g(m), D 2 m g{m) = g{m + 1) - 2g(m) + g{m - 1), m G N. (3.18) 
3.3 Bethe ansatz 

To evaluate (3.16) actually, we need to investigate the leading eigenvalue and the corre- 
sponding eigenstate. Here we derive a general formula describing the eigenvalues and the 
corresponding eigenstates through the solutions to a certain algebraic equation called the 
Bethe ansatz equation. 
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Let us define the 'vacuum state' |0) as 



|0) := 



t M/o \Utt \Utttt VO 



1 W 2 W ATi \V JVi+l \ u / JVi+2 



Ni factors 

J) «(?L~® ■■■»(! 



ATi+3 V x / ATi+4 V- 1 -/ 2iVi+2 



Ni factors 

A /o\ /o 



/ nil. I q V 1 / OAL V 1 



2ATi+3 W 2jVi+4 W yV +2iVi+2 



(3.19) 



iVo factors 

Obviously (3.19) is an eigenstate of the (twisted) quantum transfer matrix T K (X). Explicitly 
r (S (A)|0) = (a(A) + «d(A))|0), A(A)|0)=a(A)|0), £>(A)|0> = d(X)\0), (3.20) 



where 



a(A) = |ih(A^)/ 



shA 



sh(A + ei) sh(A-ei) | 2 f sh(A - e ) 1 2 sh(A - e) 



d ^ { sh(A + ei + rf) sh(A - e x + rj) ) \ sh(A - e + 77) J sh(A-e + r/) e ' T ' ^ 3 ' 21 ^ 

In the framework of the algebraic Bethe ansatz, the vector |{A} K ) constructed by the multiple 
action of B(X), namely |{A K }} = Y\!f=i B{Xj )|0), is an eigenstate of T K (X) if the rapidities 
{XjjjL-t satisfy the following Bethe ansatz equation: 

a(A*) _ ™ sh(A| - Ag + 7/) 
The corresponding eigenvalue is given by 



il sh(A - A? — 77) sh(A - A? + 77) 

a <*. -) = °« n + ^ n I (A _' A .) ■ 

.7=1 j .7=1 j 



(3.23) 



Hereafter we restrict ourselves on the case n = 1 in (3.22) and simply write the roots 
and the eigenstate corresponding to the leading eigenvalue Ao(0)(= Ao(0, 1)) as {Aj}^^ 
(N = Nq + 2N\ + 2) and |{A}), respectively. To make the analysis possible even in the 
Trotter limit N — > 00, we utilize a powerful method as in [21-23]. Let us consider the 
following auxiliary function 



N_ 
1 



1"4n m iV = iV + 2 Wl +2, (3.24) 



8 



which associates the Bethe ansatz roots {Xjjfl^ with zeros of 1 + a(A). By studying the 
analyticity properties of the auxiliary function, one sees o(A) satisfies the following nonlinear 
integral equation: 

h sh(A + 77) sh(A - e) N sh(A + rj) sh(A - e Q ) 
In a(A) = - - + In ——r + In ■ 



T sh(A) sh(A - e + rj) 2 sh(A) sh(A - e + v) 
^ N ± ^ sh(A + 77) sh(A + ei) Aq ^ sh(A + 77) sh(A - £1) 



2 sh(A) sh(A + £1 + 77) 2 sh(A) sh(A - e 1 + 77) 

dco sh(2r/)ln(l + a(w)) 
c 27ri sh(A — w + 17) sh(A — uj — rj) 



(3.25) 



Here the contour C is taken, for instance, as a rectangular contour whose edges are parallel 
to the real axis at ±7ri/2 (respectively ±77/2) and are parallel to the imaginary axis at ±77/2 
(respectively ±00) for the off-critical regime A = ch 77 > 1 (respectively for the critical regime 
< A = chrj < 1) (see figure 4 for a pictorial definition). In (3.25) the limits Nq,N± — > 00 
and £ — > can be taken analytically. We thus obtain 

1nnm _ h 2Jsh 2 (r ? ) f du S h(27y)ln(l + aH) 

{ ' T T sh(A) sh(A + rj) J c 2vri sh(A - u + rj) sh(A - u - rj) ' { ' 

which is exactly the same as that in [11]. For later convenience, here we also introduce 
another auxiliary function a(A) = l/a(A) satisfying the following nonlinear integral equation 
in the limits 7V"o, Aq — > 00 and e — > [11]: 

= h _ 2Jsh 2 (77) y d^> S h(2r/)ln(l + a(a;)) 

l j T Tsh(A)sh(A-7?) i c 2vri sh(A - w + 77) sh(A - - 77) ' l ' J 

By this auxiliary function a(A), the leading eigenvalue of the quantum transfer matrix 
related to the free energy density / by / = — TlnAo(O) is expressed as a single integral form 

lnAo( 0) = — + f ^ftfeWl + 'M) =-]L- [ ^jWlM. ( 3 . 28 ) 
w 2T 7 C 2vri sh(w) sh(u; + 7/) 2T 7 C 2vri sh(w) sh(u; - rj) y ' 

Differentiating (3.28) with respect to h, one has the magnetization (multiplied by a factor 2), 

duj sh(r])dha(uj) 1 m f du> sh(r])dha(u>) 

sh(w) sh(w — 77) (1 + a(o>)) ' 
(3.29) 



/ Z \ = i + t [ duJ sh(rj)d h a(u;) = _i_ T f 

Jc vri sh(cj) sh(w + 77) (1 + a(u;)) 7 C 7ri 



4 Multiple integral representation 

In this section, using the method developed in [11], we derive a multiple integral representa- 
tion of the longitudinal dynamical correlation function (of(0)<7^ +1 (i)) at finite temperature. 
Utilizing the relation T K (0)T K (rj) = na(rj)d(0), which can easily be obtained from (3.23), one 
expresses the operator m (3.13) as 

= K-^-Hrr +1 (^)^(-£i)]T5 L+1 (77)[T- m - 1 (£)T-^(-£ 1 )]T-^- 1 (r7). (4.1) 
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Figure 4: The integration contours for the off-critical regime A > 1 (a) and for the critical 
regime < A < 1 (b). 



To analyze the above quantity, we conveniently introduce a set of parameters £i, . . . ,£m+m 
(m = Ni/2 + m + 1, m = iVi/2 + 1) located inside C and define 

(xi, ■ ■ ■ , ^m+m) O^li • • • j -^l) • • • j -^m) 

= (£i + £, . . . ,£ m+ i + £>Cm+2 - £l, ■ ■ ■ ,£m _ e l5 Cm+1 + ^5 • • • > £m+m + 

(4.2) 

It follows that 

(tf^J*) = hm_ <M«|{£}) 

sl )"')Sm+m 

<M«l{f}) = « 3 M 

where {A} and |{A}) are, respectively, the Bethe ansatz roots and the eigenstate corresponding 
to the leading eigenvalue Ao(0) (see the previous section). Note that the dual vector ({A}| is 
constructed by the multiple action of C(A) on the state (0| which is the transposition of the 

vacuum state |0). Namely ({A}| = (0| Ufli c ( x j), where ^ = ^0 + 2^1 + 2. 

Expression (4.3) is formally similar to equation (75) in [11]. The essential difference only 
lies in the definition of the parameters {x} (4.2). In the present case, (i) the number of 
parameters {x} and elements of {x} explicitly depend on the Trotter number N±; (ii) in the 
homogeneous limit £j — > 0, the elements of {5;} converge to rj which is outside the contour C. 
On the other hand, for the static correlation function [11], the number of the parameters is 
equal to the distance of the correlator i.e. m. Furthermore all the parameters converges to 
(i.e. inside C) in the homogeneous limit £j — ► 0. 

To extend the formulation in [11] to the time-dependent case, we first introduce the 
following lemma, which is still applicable to the present case. 
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Lemma 1 [11]. $at(k|{^}) (4.3) has the following representation as a sum over partitions of 
the Bethe ansatz roots {A}, and of the inhomogeneous parameters {x} (4.2) (or equivalently 

{£})■ 

* t -™V^ V- ^ + |({A + }|{x+})^ +| ({A+}|{x}) 

$Jv( K {£}) = « > > + , n^H , (4.4) 

;s {AH r M A- } n£' a^jnJrJa + ^-)) 

{ 2; }={x+}U{x-} 
|A+|=|x+|=ti 

where |A|, |x|, etc denote the number of elements of {A}, {x}, etc. Tae too functions 
^«({A + }|{a; + }) and Z n ({A + }|{x}) are respectively defined by 



Yn({X + }\{x + }) = H 
3=1 



m+m—n 



z n {{\+}\{x})= n 



b+(A+) " sh(A+ - x+ + rj) sh(A+ - x+ - r,) 
b' + (4) M sh (4 -4 + V) sh(A+ - A+ - r?) 
detM(A+,x+)detG(A+,x+), 



l + «a(xj) J] 



where 



rn+rn 



sh(A - fi + n) 
sh(A — /i) ' 



(4.5) 



— _u a. -l u_ A sh(A+ - A+ - 7/) sh(A+ - X+ + 7/) 

and G(A, x) zs i/ie solution of a linear integral equation 

duo sh(2n) G(oj,x) 



(4.6) 



G(A,x) =t(x,X) + J 



C 27ri sh(A — to + 77) sh(A — w — 77) 1 + a{u) 





.{ 



/or x E {x} 



sh(2q) 



+ 



sh(2q) 



with 



sh(\—x+r))sh(\—x—r)) sh(A— sh(A— x+2rj) 

sh(n) 



i+al-n) } T+a(x} forx£{x} 



t(A,/z) = 



(4.7) 



(4.8) 



sh(A — /j) sh(A — fi + 77) 

Note that the range of the variable x in (4.7) is extended from inside the contour C to outside 
C by analytic continuation. 

In order to proceed to the dynamical case, let us modify (4.4) in Lemma 1. After simple 
calculations, we obtain 



<M«itf» = EE E E E 

n=0 K=0 {\}={\+}u{\-} {x}={x+}U{x-} {x}={x+}U{x-} 
\\+\=n=n+n \x + \=n \x + \=n 

_ K % + |,|J + |({A + }|{x + })% +M j + |({A+}|{x})% +us+ |({A+}|{x}) 



nl Si 1 a'(At) nPa + nHa + a(x7)) 



(4.9) 
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where n and n, respectively, denote the number of elements of {x + } and {x + } (i.e. n 
n = \x + \), n = n + n, and 



= \x+\ 



Y nfi {{\+}\{x+}) 



n 



b + (A+)b-(A+) * sh(A+ -x+ + f]) JL sh(A+ - £+ - r/) 
b' + (*+)b'_(x+) W sh(x+ - x+ + r?) 1 = 1 sh(x+ - x+ - r?) 



n 

><n 

fa— n 



sh(A+ 



77) 



sh(A+-A+-r?) 



det M(Xj,Xk) det G(Xj,Xk), 



Zn,n{{\ + }\{X}) = ]J 



1 + «0(XT)IJ 



n 



Z^({A + }|{x}) 
Here b±(A) and b±(A) denote 

fa 

b ± (A) = n 



/(*+ xj)f(xj, A+) 



(4.10) 



sh(A - x fc ) 
sh(A - x fc ± 77) ' 



ma) = n 



sh(A - x k ) 



(4.11) 



The remaining task is to replace the sums over the partitions of the set {A} and {x} in 
(4.9) with a certain set of contour integrals, where the Trotter limit will be taken analytically. 
Consider first the partitions for the set of the Bethe ansatz roots {A}. Let f(cJi, ■ ■ ■ ,co n ) be 
analytic on and inside the contour C, symmetric with respect to n variables LOj, and zero 
when any two of its variables are the same. The poles of the function 1/(1 + a(uj)) inside C 
are simple poles at u> = Xj with residues l/a'(Xj) (j = 1, . . . ,N/2). Hence the following is 
valid. 



duj-i 



2yri(l + a(uj)) 



f(ui, ...,w n ) 



E 



/(A+ ...,A+) 



{A}={A+}U{A-} n"=lO' 
|A+|=n 



(4.12) 



Note here that the relation similar to the above is also holds for a(uj). If the summand in (4.9) 
is considered to be a function of n variables {A^}" =1 , one finds it has simple poles at A^~ = x^ 
and X~j = x^ — rj. Since the parameters {£} (4.2) can be chosen arbitrary values inside C, 
we choose {x} such that the two sets of parameters {x} and {A} are distinguishable. Then 
there exists a simple closed contour surrounding the Bethe ansatz roots {A} but excluding 
{x} (see figure 5). Let C — T — T be such a contour, where T (respectively T) encircles {x} 
(respectively {x}). Applying (4.12) into (4.9), we obtain 



E 



Ffi 1 s({A + }|{x+})Z fi)S ({A+}|{x})Z fiiS ({A+}|{x}) 



{A}={A+}u{A - } uu «'(\ + ) nr=7(i + a(xj)) nr=T(i + s(2D) 

|A+|=n 



n 



(c-r-r)" =1 



da;, 



2vri(l + 0(0^)) 



%({ w }|{x+})%({ w }|W)z fiiS (M|{x}) 



nr=T(i+^-))nr=T(i+^-)) 



(4.13) 
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Figure 5: The integration contour C — T — T corresponding to the off-critical regime A > 1. 



Next step is to reduce the integrals along the contour C— T— T to those along the canonical 
contour C. We first consider the integrals on T. Because the integrand in (4.13) is symmetric 
with respect to {u;}, we can divide the integral 



(c-r-r)« , =1 



2vri 



£(-i)< 



k=l 



TT — / IT 

(c-r)»-* fJi 27ri M p x 



n—k+j 



2vri 



(4.14) 



where n = n + n. Note that the sum over k is restricted to k < n, since \x\ = n and the 
integrand vanishes when any two of {uj} is the same. Noting that, inside T, G(ujj,Xk) (4.7) 
has simple poles Uj = with residues —1, and the poles at Uj = Xk for M(uj,x^) (4.6) are 
canceled by simple zeros of b+(w), we have 



Li 



i=i 
fc! 



n—k+j 



2ni(l + a(uj)) 



yji,n(H}p =a |{x + })^ )S ({ Wi }p =1 |{x})^, s ({ Wi }7 =1 |{x}) 



^ y a - fc ,n(H-}j=ll{^ + » JJ 

{x+-}U{x++ }={£+} n j =l( 1 + a ( X j h )) j=l 
{x++}U{£+}={x++} 
|2+-]=fe 



n— fc 



l + zc-'sfij-jj] 



><n 



n— A; 



1 + Kd(xj) Y\ 



f{x j ,x l )f(ui,x j ) 



\ f{xJ + ,X j )f(x j ,U|) 



n 

i=i 



f( x t + > x j )f( x j > w i) 
i f{ x J> x i' + )f( u; l> x j) 



^ fc /(£}-, x++)/(^,x+-) 

K fJM + ^-)/(2+-,^)_ 



(4.15) 
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Performing a resummation as in [11], one obtains 



<m«ik»=EE E 



E 



_to— m— n 



71 



1 



ii=os=o{j+}u{j-}={x} {£+}u{£-}={£} j'=i -n ^ 

|3; + |=n |x+|=n 



= / - n 

'(c-r)» ^ 



do;jb_(a;j)b+(a;j) -A- sh(u;j — x k — rj) -A sh(wj — x\ + 77) 



n 



n 



2vri(l + a(Wj)) sh(x+ -x\ -r\) ^ sh(x+ - x+ + 77) 
x det M(uij,x^) det G(ujj,x k ) 



R 1 + a(5£ ) 2 



1 + Ka(x k ) J] 



f(uj>x k )f(x k ,xp 

1 f@k>Uj)f( X j'Xk) 



(4.16) 



where 



M( Wj ,x+)=i(x+ a;,) J]- 



sh(ajj — X^ — 7]) 



L sh(u^ -ui-rj) 



sh(ujj — X^ + 7]) 



s h(wj - + 77) 



(4.17) 



Now we would like to consider the integrals along the contour T. Utilizing the relation 
similar to (4.14), we can separate the C-integrals from the T-integrals. The unwanted terms 
including a(xj) etc, however, cannot be eliminated by actually performing the T-integrals, 
since the poles inside F are not {x} but {x — 77}. Nevertheless, we observe that all the 
unwanted terms vanish in the homogeneous limit £j — > 0, because the auxiliary function o(A) 
has poles (respectively zeros) of order Nq/2 + N\ + 1 at A = (respectively A = 77). From 
the relation lim/^j^/o} G(ujj,Xk) = — lim/^i^o} G(ujj,Xk — rj) which is derived by adapting 
a(0) = 00, a(rj) = and x) = t(x — 77, w) to (4.7), we obtain 



^(»|(°» = {{) %EE E 



E 



1 



n! - LJ - b' fr+W (V +N l 

n=0n=0 {x+}U{x-}={x} {x+}U{x~}={x} 3=1 ~ V j ^+^.7 ^ 



in 



|a; + |=n 

dLOjb^(ujj)b + (ujj) -A sh(wj — x^ — 77) 
2vri(l + a(cjj)) ^ sh(x+ - x£ - 77) 

Here the function R(cjj, x£) is defined by 

G(u)j,xl) 



detM(w j ,x+)deti?(u; i ,x+). (4.18) 



R(ujj,xJ 



-KG{uj,x+ -ri) nr=i 



sh(a:+ -a^-77) sh(a:+ -3;+ +7?) 



for x^" G {x} 
for x^ € {x} ' 



(4.19) 



and we have used fh — m = m and n + n = n. The integrand of (4.18) is a symmetric 
function of {x + } and vanishes at x+ = x^ and x+ = xjj" . Thanks to this together with the 

fact that l/b_(A) (respectively l/b + (A)) has simple poles at A = x (respectively A = x), we 
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can directly apply (4.12) to (4.18). Thus, we finally arrive at 

do 



m+m m — n n 

$ n (k\{0}) = lim V ^tTI 
sh(uj -Ck~v) 



r{o}ur{,} 2vrib_(Ci)b+(Ci) Jc 27ri (l + a(wj)) 
detM(a; J -,a)deti?(a; j ,a), (4.20) 



where the contour T{0} U T{n} surrounds the points and rj and does not contain any other 
singularities. In (4.20), the Trotter limit Nq — > oo, iVi — ► oo and the limit e — > can be 
taken analytically. In this limit, the auxiliary function a(A) is the solution to the nonlinear 
integral equation (3.26). Substituting (4.2) and (3.3), we easily obtain that 



lim lim b_(A)b+(A) = lim 



sh(A + /?i/JVi) 


2 


sh(A — rj) 


2 


shA 


sh(X + p 1 /N 1 -r ] )_ 




sh(A) 




sh(A-??) 



exp 



— ite — 2) — ( 



(4.21) 



where e(A) and p{\) are related to the bare one-particle energy and momentum, respectively. 
They are explicitly given by 



e(A) 



2Jsh 2 ?y 



8 h(A+3)sh(A-2)' ""Vsh(A-f)r 
From (4.20), (4.3) and (3.16), we end up with the following theorem. 

Theorem 1 The generating function of the longitudinal dynamical correlation function has 
the following multiple integral representation: 



(n!) 2 

n=0 V ; j=l 



dC, 



duJn 



r{o}ur{r,} 2vri 27ri(l + a(wj)) 



t~t sh(u; i -Ck-V) 



x J^e"^-^-^-^^^^^ (4.23) 

j'=i 

in £/ms expression the function a(A) is i/te solution to the nonlinear integral equation (3.26); 
e(A) and p(A), respectively, denote the bare one-particle energy and momentum (4.22); the 
function M(u,Q and R(lo,Q are, respectively, defined by 



M(„c)^(^)n :£:!:il + ^ong 



sh(a; -Cl+V) 



J= i sh(w -ui-rj) 



KU[U,(, V)lll=l sh(C-^+r,)sh(C-0-r?) /0r ^ 



a; - uji + n) ' 

for C ~ 
~ 77 ' 



(4.24) 
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where t(X,fi) = sh(??)/(sh(A — fj,) sh(A — n+rjj) and G(X,Q is the solution of the linear integral 
equation 

G(\, C) = t(C, A) + f ^ --^V (4.25) 

v ' s; vs ' ; Jc 2th sh(A - u + ri) sh(A - u - rj) 1 + o(w) V ' 

The contour C is the canonical contour (see figure 4) and T{0} U T{r]} encircles the points 
and rj. 

Thus the time-dependent longitudinal correlation function (ci(0)cr^ +1 (i)) can be evalu- 
ated by inserting the generating function (4.23) into (3.17), and using the magnetization 

which is derived from (3.26), (3.27), (3.29) and (4.25) [11]. 

5 Special cases 

Here we comment on some special cases derived from the multiple integral representation 
(4.23) for the dynamical correlation function (erf (0)cr^ +1 (t)). 

5.1 Static limit t = 

First let us consider the static limit t = 0. Due to the factors e((j — rj/2), the integrand in 
(4.23) has essential singularities located on Q = and 77. In the static limit these essential 
singularities vanish, and hence the integrals along the contour T{r/} disappear. Because the 
remaining part of the integrand has poles of order m at the points Cj = 0) the integrals on 
the contour T{0} vanish for n > m. Namely the sum over n is restricted to n < m. The 
resultant expression reads 



a.K°)=E7^2 n 

n=0 y '> j=l 



d(j ( sh(Q — rj) \ m f dujj ( sh(uj 



r{o} 2vri [ sh(Cj) J Jc 27ri(l + a(uj)) { sh(uj - rj) 



x II 6/ ? i detM(^,C fc )detG(^,C fc ), (5.1) 



agreeing with that in [11]. 



5.2 Zero-temperature limit T = 

For direct comparison with the previous result for T = [9], it is convenient to intro- 
duce another expression of the generating function. Considering the correlation function 
*) cr i(0)) which is equivalent to (o~i(0)a^ l+1 (t)), and utilizing the technique described 
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in the previous section 3 , one may obtain 



n 



(n!) 2 

n=0 v ; j=l 



duj-i 



r{_r,}ur{o} 2ir\ Jc 27ri (! + <*(^)) 



TT M^j -Ck + v) 

^ MC 3 -Ck + v) 



x Y[ c i *( e (0+I)-^K'+I))+™(p(C 3 +!)-PK+!)) det M(uj,Ck) det 

where Jj(A) and C) are defined by 

'sh(A-gV 



(5.2) 



p(A) = iln 



sh(A+2) I ' 



^(w,C) = S ..-1 



for C ~ 



-K^Cld,! f + n\T\ n sh «-^l+V) M<-Cl-V) for C rsj 77 

^^'S + 7 /Vll/=l B h(C-u;j-7;)sh(C-Ci+r/) ^ ^ 



(5.3) 



Note that G(A, £) is the solution of the integral equation (4.7) which is also written in terms 
of a(A) if A and £ are located inside C: 



G(\x) = -t(\X) - j c 



sh(2r/) 



G(u,C) 



C 2m sh(A — lo + rj) sh(A — uj — rj) 1 + a(w) ' 



(5.4) 



Here we restrict ourselves on the off-critical case A > and set r/ < as in [9]. Note that 
we can also treat the critical case |A| < 1 by just changing the definition of the integration 
contour as in figure 4. 

Shifting the variables in (5.2) by ujj — > Uj—v/2 and Q — ► Q —rj/2, we consider the integrals 
on the contour T{-r]/2} U T{r]/2} and -C UCh. Here C and Ca denote C = [— vri/2, vri/2] 
and Ca = [r//2 — 7ri/2,r?/2 + vri/2], respectively. By close analysis of the auxiliary function 
o(A) for h > and 77 < at the zero-temperature limit T — > 0, one finds 



1 t-o JO for A € C U Cr \ C 
l + a(A) " I 1 for A e £ 



rj 77 
-9h + -,% + - 



(5.5) 



where the 'Fermi point' is an imaginary number depending on h. Substituting this into 
(5.4) and appropriately shifting the variables, we have 

° \ X ~ 2 ' C " 2 ) ~ " t(A ' + 7_ £ 2^sHX-co + V ) S h(X-co- V Y (5 - 6) 

Comparing above with equation (2.16) in [7], one identifies G(X, Q as the density function 
P(A,C)= 

G(A-|,C-|) =2vrip(A,C). (5.7) 



In this case, one needs to modify the definition of the spectral parameters in the quantum transfer matrix 
as A ± Ej — ► A and —A — » —A ± Ej , where Sj = e, £0, £i- 
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Inserting (5.7) and (5.5) into (5.2), we arrive at the expression for the generating function at 
T = 0: 



h{±r ] /2} 2vri J_ c 



sh (o -Ck+v) 



n 



x ei^efcO-eK^+imbfcO-pK)) det M (uj,Ck) det n(ujXk) 

3=1 



with 



/»(w,C) 

The above expression reproduces equation (6.17) in [9]. 



for C ~ 
^RC^PIOTI tor ^ _r? 



(5.8) 



(5.9) 



5.3 XY (free fermion) limit A = 

Along the method described in [9] (see also [24,25] for the static case), we would like to 
study the XY limit, where (cf (0)cr^ +1 (t)) can be written as a product of single integrals. Set 
rj = vri/2. Then the kernels of the integral equation in (3.26) and in (4.25) are equal to zero. 
Hence 



a(A) = exp 



h 



1 4iJ 



T Tsh(2A)J' ~ v "'^ sh(2(A-C))' 
Shifting the variables Q — > Q + vri/4 and ujj — > Uj + 7ri/4 in (4.23), one easily sees 



(5.10) 



n=0 v ' 7 j=l L 



Jr{±7ri/4} 27ri J a 



27ri(l + o(wj + f i)) 



x n e i *( e ^')- e ^)) +im ( p ( c: ^- p ^)) det M^Xk) det R^j^k), (5.11) 
i=i 

where C = —Co U C_ 7ri / 2 ; Co = [ — oo, oo]; C_ 7r i/ 2 = [—vri/2 — oo, — 7ri/2 + oo]. In this case the 
function R(u, C) and M(w, C) are reduced to 



R(u,0 



Bh(2(w-C)) 
2 



for £ ~ 7ri/4 



2(/e- 1) 



for C ~ -7ri/4 ' ® ' s1i(2(kj - ()) ^ ch(u - wj) 



n 



ch(w - 0) 



(5.12) 



k '*Bh(2( W -0) 

In the above expression, we notice that M(uj, Q is factorized and has the factor k — 1, which 
significantly simplifies the integral representation. After taking the second derivative with 
respect to k and setting k = 1, one observes all the terms n > 2 vanish. Firstly we consider 
the case n = 2. Substituting (5.12) into (5.11), we extract the term corresponding to n = 2 

(2) 

(written as Q K (m,t)). After differentiating with respect to k and setting k = 1, one has 

2 



\k=1 



32tt / 



n 



Jc 1 + a(uj + f ) 7r{±7ri/4} 



det 



e We(Cfc)+imp(Cfe) 

sh(o; j - Ck) 



det 



g — it e (oj j ) — imp (w ) 

sh(o; j - ( k ) 



(5.13) 
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The integral on r{±7ri/4} can be evaluated by considering the residues outside the contour 
r{±7ri/4} i.e. at the points £j = uj k . This leads to 



a 1 + a(uj + f ) 



det 



1 _ e it(e(wj)-e(w k ))+im(p(wj)-p(w k )) 



sh(uJj - UJ k ) 



Changing the variables cosh(2u;) = — 1/ cosp and identifying 
ch(2w) 



duj = dp, e{to) = 4 J cos p, 



y/dh(2LJj)dh(2u) k ) 



one obtains 



-1 



«=1 167T 2 



n 



/*7T 




/ dpji9(pj) 


det 







s\i{uj - u k ) smM Pj -p k ) 



I — e 4itJ(cosp J -cosp fc )+im(p J '-p fc ) 

sin ±(pj -p k ) 



(5.14) 
(5.15) 



(5.16) 



where $(p) is the Fermi distribution function 

0(p) = 



1 



l + exp[-^ + ^£] 



(5.17) 



Taking the lattice derivative over m, we obtain the following relation, without explicit eval- 
uation of the multiple integral: 



2D 2 m d 2 K QW(m,t) 



1 



K=l 



r d P $(p) 

J —IT 



1 



7T- 



dpd(p)e 



4itJcos p+imp 



(5.18) 



Next let us compute the term Q^\m,t) corresponding to n = 1. Its explicit form is 



e«(»M) = 



K" 1 " 1 ^- 1) 

4^2 



dC + k 



X 



r{-7ri/4} 
Ic l + a 



do; 



r{7ri/4} 
e it(e(C)-e(o;))+im(p(C)-p(a;)) 



/C i + a(a; + f ) 
Evaluating the integral on r{±7ri/4}, one has 

e it(e(C)-e(oj))+im{p(0-p(oj)) 



sh 2 (u-C) 



(5.19) 



/ dC+ / dC 



It immediately follows that 
2DldlQW(m,t) 



sh 2 ( W -C) 



2ir[te'(u) +mp'(u)]. (5.20) 



da; 



jt(e(C)-e(oj))+im(p(0-p(u>)) 



k=i vr 2 7r{-7ri/4} Jc> 1 + a(u + f ) ch(2w) ch(2C) 



4 f duj ,. , 



(5.21) 
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The contour T{— 7ri/4} in the first term can be replaced by C. After changing variables as 
in (5.15), one arrives at 



dp$(p). 



ZD m dlQ£Xm,t) =X T dp$(p)e- 4itJcosp -' imp dqe 4itJcosq+imq - - T 

K=l 7T 2 J_ 7T TT J_ n 

(5.22) 

Sum up (5.18), (5.22) and 2Dl l d K Q ( -°\m, t)\ K=1 = 4 which is trivially obtained from Q^(m, t) 
K m . Then inserting the result into (3.17) and using the magnetization (see (4.26)): 

{a z } = -1 + 1/ dp#(p), (5.23) 



we finally obtain 

1 



K(o)< +1 W)=K) 2 -- 2 



r dpi9(p)e 4it 

J —TT 



4it J cos p+\mp 

+ \ dp$(p)e- 4itJcosp -' imp r dqe AitJcosq+imq . (5.24) 

The above expression coincides with the known result as in [15] 4 . Of course, (5.24) can also 
be derived by starting from the generating function defined in (5.2). 

6 Conclusion 

In this paper, we have derived a multiple integral representation for the time-dependent 
longitudinal correlation function (af (0)a^ t+1 (t)) at any finite temperature and finite magnetic 
field. The formula reproduces the known results in the following three limits: (i) static limit, 
(ii) low-temperature limit and (iii) XY limit. 

It will be very interesting to consider the following problems which still remain open. First, 
it is well-known that the long-distance asymptotics of correlation functions at the low-energy 
region (namely T = or T <C J) in the critical regime can be derived by a field theoretical 
argument (see [26] for example). In contrast, our multiple integral representation (4.23) is 
valid for any finite temperature and interaction strength. The exact computations of the 
asymptotic behavior from (4.23) beyond the field theoretical predictions are of importance. 

Second, how do we evaluate and extract the long-time asymptotics of the correlation 
function (af (0)a^j, 1 (i)) at infinite temperature? In this case the auxiliary function a(A) in 
(4.23) becomes quite simple: o(A) = 1. This problem is of interest in relation to the issue of 
spin-diffusion in the spin-1/2 Heisenberg XXZ chain (see, for example, [27-30]). 

The third is how to apply our formula to the calculation of crucial physical quantities 
such as the dynamical spin structure factor [31-36], which can be actually measured by a 
neutron scattering experiment. In relation to this problem, finally we would like to comment 
on the form factor expansion. In fact, the multiple integral representation for T = (see 
(5.8)) is directly connected to the form factor expansion [8,9]. On the other hand, for 
the finite temperature case, it is not clear whether (4.23) has a connection with the form 



4 In fact the correlation function {o'm+i(t) cr i(0)) is considered in [15]. The result is the same as that of 

W(o)< +1 (-t)). 
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factor expansion, since (4.23) is derived by the quantum transfer matrix acting not on the 
quantum space but on the auxiliary space. The form factor expansion is an important tool 
to investigate the dynamical properties of the system, and therefore explicit expressions at 
finite temperature are also desired. 

Acknowledgments 

The author would like to thank J. Sato and M. Shiroishi for fruitful discussions. He also 
acknowledges with thanks the hospitality of the theory group of the Australian National 
University, where part of this work was done. This work is partially supported by Grants- 
in-Aid for Young Scientists (B) No. 17740248, Scientific Research (B) No. 18340112 and (C) 
No. 18540341 from the Ministry of Education, Culture, Sports, Science and Technology of 
Japan. 

References 

[I] M. Jimbo, K. Miki, T. Miwa, and A. Nakayashiki, Phys. Lett. A 168 (1992) 256. 

[2] M. Jimbo and T. Miwa, Algebraic Analysis of Solvable Lattice Models, (American Math- 
ematical Society, Providence, RI, 1995). 

[3] A. Nakayashiki, Int. J. Mod. Phys. A 9 (1994) 5673. 

[4] V.E. Korepin, A.G. Izergin, F.H.L. Essler, and D.B. Uglov, Phys. Lett. A 190 (1994) 182. 

[5] M. Jimbo and T. Miwa, J. Phys. A 29 (1996) 2923. 

[6] N. Kitanine, J. M. Maillet and V. Terras, Nucl. Phys. B 567 (2000) 554. 

[7] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, Nucl. Phys. B 641 (2002) 487. 

[8] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, Nucl. Phys. B 712 (2005) 600. 

[9] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, Nucl. Phys. B 729 (2005) 558. 

[10] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, hep-th/0505006. 

[II] F. Gohmann, A. Klumper and A. Seel, J. Phys. A 37 (2004) 7625. 
[12] F. Gohmann, A. Klumper and A. Seel, J. Phys. A 38 (2005) 1833. 

[13] F. Gohmann and N. P. Hasenclever and A. Seel, J. Stat. Mech. 0510 (2005) P10015. 
[14] H. E. Boos, F. G6hmann, A. Klumper and J. Suzuki, J. Stat. Mech. 0604 (2006) P04001. 
[15] F. Colomo, A.G. Izergin, V. E. Korepin, V. Tognetti, Theor. Math. Phys. 94 (1993) 11. 
[16] N. Kitanine, J. M. Maillet and V. Terras, Nucl. Phys. B 554 (1999) 647. 
[17] J. M. Maillet and V. Terras, Nucl. Phys. B 575 (2000) 627. 



21 



[18] F. Gohmann and V.E. Korepin, J. Phys. A 33 (2000) 1199. 

[19] M. Suzuki, Phys. Rev. B 31 (1985) 2957. 

[20] M. Suzuki and M. Inoue, Prog. Theor. Phys. 78 (1987) 787. 

[21] A. Klihnper, Ann. Phys., Lpz. 1 (1992) 540. 

[22] A. Klumper, Z. Phys. B 91 (1993) 507. 

[23] C. Destri and de H. J. de Vega, Phys. Rev. Lett. 69 (1992) 2313. 

[24] N. Kitanine, J. M. Maillet, N. A. Slavnov and V. Terras, Nucl. Phys. B 642 (2002) 433. 

[25] F. Gohmann and A. Seel, Theor. and Math. Phys. 146 (2006) 119. 

[26] V.E. Korepin, N.M. Bogoliubov and A.G. Izergin, Quantum Inverse Scattering Method 
and Correlation Functions (Cambridge University Press, Cambridge, 1993). 

[27] K. Fabricius, U Low and J. Stolze, Phys. Rev. B 55 (1997) 5833. 

[28] K. Fabricius and B.M. McCoy, Phys. Rev. B 57 (1998) 8340. 

[29] S. Glocke, A. Klumper and J. Sirker, cond-mat/0610689. 

[30] J. Sirker, Phys. Rev. B 73 (2006) 224424. 

[31] D. Biegel, M. Karbach and G. Miiller, Europhys. Lett. 59 (2002) 882. 

[32] J. Sato and M. Shiroishi, J. Phys. Soc. Jpn. 73 (2004) 3008. 

[33] J.-S. Caux, R. Hagemans, J.-M. Maillet, J. Stat. Mech. (2005) P09003. 

[34] J.-S. Caux and J. M. Maillet, Phys. Rev. Lett. 95 (2005) 077201. 

[35] R. G. Pereira, J. Sirker, J.-S. Caux, R. Hagemans, J. M. Maillet, S. R. White, I. Affleck, 
Phys. Rev. Lett. 96 (2006) 257202. 

[36] J.-S. Caux and R. Hagemans, J. Stat. Mech. 0612 (2006) P12013. 



22 



